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ABSTRACT 

Context. Trans-Neptunian objects (TNO) represent the leftovers of the formation of the Solar System. Their physical properties 
provide constraints to the models of formation and evolution of the various dynamical classes of objects in the outer Solar System. 
Aims. Based on a sample of 19 classical TNOs we determine radiometric sizes, geometric albedos and beaming parameters. Our 
sample is composed of both dynamically hot and cold classicals. We study the correlations of diameter and albedo of these two 
subsamples with each other and with orbital parameters, spectral slopes and colors. 

Methods. We have done three-band photometric observations with Herschel/PACS and we use a consistent method for data reduction 
and aperture photometry of this sample to obtain monochromatic flux densities at 70.0, 100.0 and 160.0 yum. Additionally, we use 
SpitzerfMlPS flux densities at 23.68 and 71.42 yum when available, and we present new Spitzer flux densities of eight targets. We 
derive diameters and albedos with the near-Earth asteroid thermal model (NEATM). As auxiliary data we use reexamined absolute 
visual magnitudes from the literature and data bases, part of which have been obtained by ground based programs in support of our 
Herschel key program. 

Results. We have determined for the first time radiometric sizes and albedos of eight classical TNOs, and refined previous size and 
albedo estimates or limits of 1 1 other classicals. The new size estimates of 2002 MS4 and 120347 Salacia indicate that they are among 
the 10 largest TNOs known. Our new results confirm the recent findings that there are very diverse albedos among the classical TNOs 
and that cold classicals possess a high average albedo (0.17 ± 0.04). Diameters of classical TNOs strongly correlate with orbital 
inclination in our sample. We also deteiTnine the bulk densities of six binary TNOs. 

Key words. Kuiper belt - Infrared: planetary systems - Techniques: photometric 



1. Introduction 

The physical properties of small Solar System bodies offer constraints on theories of the formation and evolution of the planets. 
Trans-Neptunian objects (TNO), also known as Kuiper Belt objects (KBO), represent the leftovers from the formation period of the 
outer Solar System (IMorbidelli et al. 200 8i. and they are analogues to the parent bodies of dust in debris disks around other stars 
( |Wyatt 2008|IMoro-Martm et al. 2()08l and references therein). 

* Herschel is an ESA space observatory with science instruments provided by European-led Principal Investigator consortia and with important 
participation from NASA. 
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In addition to Pluto more than 1400 TNOs have been discovered since the first Kuiper belt object in 1992 jJewitt and Luu 19^3T l. 
and the current discovery rate is 10 to 40 new TNOs/year. The dynamical classification is based on the current short-term dynamics. 
We use the classification of Gladman ( jGladman et al. (2008) 10 Myr time-scale): classical TNOs are those non-resonant TNOs 
which do not belong to any other TNO class. The classical TNOs are further divided into the main classical belt, the inner belt 
(a < 39.4 AU) and the outer belt (a > 48.4 AU). The eccentricity limit for classicals is e < 0.24, beyond which targets are classified as 
detached or scattered objects. The classification scheme of the Deep Eplictic Survey Team (DES. IETliot et al. 1 2005) differs from the 
Gladman system in terms of the boundary of classical and scattered objects, which do not show a clear demarcation in their orbital 
parameters. Some of the classicals in the Gladman system are scattered-near or scattered-extended in the DES system. Another divi- 
sion is made in the inclination/eccentricity space. Although there is no dynamical separation, there seems to be two distinct but partly 
overlapping inclination distributions with the low-i "cold" classicals, limited to the main classical belt, showing different average 
albedo ( [Grundy et al. 2005 IBrucker et al. 1 2009). color ( [Trujillo and Brown 2 002), luminosity function JFraser et al.l 2010). and fre- 
quency of binary systems (Noll et al. 2008) than the high-i "hot" classicals, which has a wider inclination distribution. Furthermore, 
models based on recent surveys suggest that there is considerable sub-structure within the main classical belt (Petit et al. 201 I j . To 
explain these differences more quantitative data on physical size and surface composition are needed. 

The physical characterization of TNOs has been limited by their large distance and relatively small sizes. Accurate albedos 
help to correctly interpret spectra and are needed to find correlations in the albedo-size-color-orbital parameters space that trace 
dynamical and collisional history. The determination of the size frequency distribution (SFD) of TNOs provides one constraint to 
formation models and gives the total mass. The SFD of large bodies is dominated by accretion processes and they hold informa- 
tion about the angular momentum of the pre-solar nebula whereas bodies smaller than 50 to 100 km are the result of collisional 
evolution dPetit et al 2008). The SFD can be estimated via the luminosity function (LF), but this size distribution also depends on 
assumptions made about surface properties such as albedo. Consequently, ambiguities in the size distributions derived from the LFs 
of various dynamical classes are one significant reason why there is a wide uncertainty in the total TNO mass estimate ranging from 
0.01 Mearth (IBemstein et al. 2004) to 0.2 MEanh ( [Chiang et al. 1999 |l. Among the formation models of our Solar System the "Nice' 



family of models have been successful in explaining the orbits of planets and the formation of the Kuiper belt ( Tsiganis et al. 2005 
ILevison et al. 2008'). although they have difficulties in explaining some of the details of the cold and hot distributions and the origin 
of the two sub-populations (e.g. Eraser et aL 2010. Petit et al. 2011l[Batygin et al. 201 ip . 

Only a few largest TNOs have optical size estimates based on direct imaging and assumptions about the Umb darkening function 
(e.g. Ouaoar. Eras er and Brownl 2010). The combination of optical and thermal infrared observations gives both sizes and geometric 
albedos, but requires thermal modeling. Earlier results from Spitzer and Herschel have shown the usefulness of this method (e.g. 
Stansberry et al. 2008. Miil ler et al.l 2010) and significantly changed the size and albedo estimates of several TNOs compared to 



those obtained by using an assumed albedo. 

In this work we present new radiometric diameters and geometric albedos for 19 classical TNOs. Half of them have no previously 
published observations in the wavelength regime used in this work. Those which have been observed before by Spitzer now have 
more complete sampling of their SEDs close to the thermal peak. The new estimates of the 19 targets are based on observations 
performed with the ESA Herschel Space Observatory dPilbratt et al. 201 Oi l and its Photodetector Array Camera and Spectrometer 
(PACS; Pogli tsch et al. 2010 1). Other Herschel results for TNOs have been presented bv M iiller etaL (2010). Lellouch et aT] (2010) 
and Lim et alj(2010). New estimates of 18 Plutinos are presented in Mommert et al. (2012) and of 15 scattered disc and detached 
objects in lSantos-Sanz et an (2012). 

This paper is organized in the following way. We describe our target sample in Section lTTl Herschel observations in Section l272l 
and Herschel data reduction in Section l273] New or re-analyzed flux densities from Spitzer are presented in Section l24l As auxiliary 
data we use absolute V-band magnitudes (Section l23T l. which we have adopted from other works or data bases taking into account 
the factors relevant to their uncertainty estimates. Thermal modeling is described in Section|3]and the results for individual targets 
in SectionlH In Section|5]we discuss sample properties of our sample and of all classicals with radiometric diameters and geometric 
albedos as well as correlations (Section lSTTT i and the bulk densities of binaries (Section l5T2l i. Finally, the conclusions are in Section|6] 



2. Observations and data reduction 

Our sample of 19 TNOs has been observed as part of the Herschel key program "TNOs are Cool" ( Miiller et al. 20091 mainly 
between February and November 2010 by the photometry sub-instrument of PACS in the wavelength range 60-210 yum. 



2.1. Target sample 

The target sample consists of both dynamically cold and hot classicals (Table [TJ. We use a cut-off limit of / = 4.5° in illustrat- 
ing the two subsamples. Another typical value used in the literature is ; - 5°. The inclination limit is lower for large objects 
dPetit etal. 20111 ) which have a higher probability of belonging to the hot population. Targets 119951 (2002 KXu), 120181 (2003 
UR292) and 78799 (2002 XW93) are in the inner classical belt and are therefore considered to belong to the low-inclination tail of 
the hot population. The latter two would be Centaurs in the DES system and all targets in Table [T] with / > 15° would belong to the 
scattered-extended class of DES. 

The median absolute V-magnitudes (i/y, see Section l23] ) of our sample are 6. 1 mag for the cold sub-sample and 5.3 mag for the 
hot one. Levison and Stem (2001) found that bright classicals have systematically higher inclinations than fainter ones. This trend 
is seen among our targets (see Section l5.1.3b . Another known population characteristics is the lack of a clear color demarcation line 
at / a! 5° jPeixinho et al. 20081 ). which is absent also from our sample of 10 targets with known colors. 
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Table 1. Target sample. Semimajor axis a, perihelion distance q, inclination /, eccentricity e, color taxonomy, spectral slope, and 
the average absolute visual magnitudes (V or R-band) of 19 TNOs ordered according to increasing inclination. For i/y used in our 
analysis, see Table |5] 



Target 


KBO 


a" 


q 


i" 


e" 


Color " 


Spectral slope 






location 


(AU) 


(AU) 


n 






(% / 100 nm) 


(mag) 


119951 (2002 KXm) 


inner 


38.9 


37.1 


0.4 


0.05 


RR-IR' 


27.1 ± 1.0'' 


4.862 ± 0.038'' 


(2001 XR254)* 


main 


43.0 


41.7 


1.2 


0.03 




10 ± 3 '' 


6.030 ± 0.017' 


275809 (2001 QY297)* 


main 


44.0 


40.4 


1.5 


0.08 


BR 


24 ± 8-^'^'''' 


6.09 ± 0.03^ 


(2001 RZ143)* 


main 


44.4 


41.3 


2.1 


0.07 




13 ± 6''-'-^ 


6.69 ±0.10' 


(2002 GV31) 


main 


43.9 


40.0 


2.2 


0.09 






= 5.5' ±0.4 


79360 Sila (1997 CSjs)* 


main 


43.9 


43.4 


2.2 


0.01 


RR 


27.0 ± 3.0'' 


5.59 ± 0.06^"" 


88611 Teharonhiawako (2001 QTjgy)* 


main 


44.2 


43.2 


2.6 


0.02 




1±2/ 


5.97 ± 0.03^ 


120181 (2003 UR292) 


inner 


32.6 


26.8 


2.7 


0.18 




28 ±5* 


//r = 6.7' ±0.3 


(2005 EF298) 


main 


43.9 


40.1 


2.9 


0.09 


RR 




= 5.8' ±0.3 


138537 (2000 OKgy) 


main 


46.8 


40.0 


4.9 


0.14 


RR 


20 ± 3"-'-'-'"-" 


6.47 ± 0.09' "' 


148780 Altjira (2001 UQig)* 


main 


44.5 


41.8 


5.2 


0.06 


RR 


35 ± 6^-«-''J 


6.47 ± 0.1 3' s 


(2002 KW14) 


main 


46.5 


37.3 


9.8 


0.20 






5.88 ± 0.05'' 


(2001 KA77) 


main 


47.3 


42.8 


11.9 


0.10 


RR 


38 ± 3^-''-'-"'-° 


5.64 ± 0.08« 


19521 Chaos (1998 WH24) 


main 


46.0 


41.1 


12.0 


0.11 


IR 


23 + 2'^-}-"'-"-p-i 


4.97 ± 0.05«'" 


78799 (2002 XW93) 


inner 


37.6 


28.3 


14.3 


0.25 






//r = 4.8' ±0.6 


(2002 MS4) 


main 


41.7 


35.6 


17.7 


0.15 






Hr = 3.5' ±0.4 


145452 (2005 RN43) 


main 


41.8 


40.6 


19.2 


0.03 


RR-IR'^ 


23.0 ± i.r 


3.89 ± 0.05" 


90568 (2004 GV,) 


main 


41.8 


38.7 


22.0 


0.07 


BR'^' 


15 ± 3* 


4.25 ± 0.04"' 


120347 Salacia (2004 SBgo)* 


main 


42.2 


37.9 


23.9 


0.10 




12.6 + 2.0" 


4.26 ± 0.02^ 



Notes. The horizontal line marks the limit of dynamically cold and hot classicals according to our dynamical analy- 
sis using the Gladman system. * denotes a known binary system iNoU et al. 2008b . References. lAU Minor Planet 
Center, http://www.minorplanetcenter.net/iau/Ephemerides/Distant7] accessed July 2011. Taxonomic class from 

|Fulc hignoni"et al. (2008) unless otherwise indicated. <''' IBarucci et all (2011). 'Pema et al. (2010)' [Forn asier et al. (2009)| 
"^jBenecchi et al. (2009) '^' |Doressoundiram et al. (2005)| Calculated using the technique of Hainaut and Delsanti (2002). Santos-Sa nz et al.l 
(2009). Values from the Minor Bodies in the Outer Solar System database, http: //www. eso . org/~ohainaut/MBOSS Stephen Tegler, priv. 
co mm. <" B enecchi et al. (2011). |Do ressoundiram et al. (20"0 2)] <"> [D elsanti et a l. (2001)| Peixinho et al. (2004) i''' |Barucci et al. (2000) | 
[Tegler and Roma nishin (2000) Rabinow itz et al. (2007) T~^" Grundy et al. (2009) lAU Minor Planet Center / List of Transneptunian 
Objects at http : //www. minorplanetcente r .net/iau/lists/TNOs .html accessed June 2011. See Section |231 for conversion to /fy 
'"' [Romanishin and Tegler] (2005). Pema et al., in prep. '"'' [DeMeo et al. (2009)[ 



2.2. Herschel observations 

PACS is an imaging dual band photometer with a rectangular field of view of 1.75' x 3.5' with full sampling of the 3.5 m-telescope's 
point spread function (PSF). The two detectors are bolometer arrays, the short-wavelength one has 64 x 32 pixels and the long- 
wavelength one 32 x 16 pixels. In addition, the short-wavelength array has a filter wheel to select between two bands: 60-85 //m or 
85-125 /im, whereas the long- wavelength band is 125-210 jim. In the PACS photometric system these bands have been assigned 
the reference wavelengths 70.0 fim, 100.0 yum and 160.0 /vm and they have the names "blue", "green" and "red". Both bolometers 
are read-out at 40 Hz continuously and binned by a factor of four on-board. 

We specified the PACS observation requests (AOR) using the scan-map Astronomical Observation Template (AOT) in HSpot, 
a tool provided by the Herschel Science Ground Segment Consortium. The scan-map mode was selected due to its better overall 
performance compared to the point-source mode ( Miiller et al.l 2010). In this mode the pointing of the telescope is slewed at a 
constant speed over parallel lines, or "legs". We used 10 scan legs in each AOR, separated by 4". The length of each leg was 3.0', 
except for Altjira where it was 2.5', and the slewing speed was 20"^"'. Each one of these maps was repeated from two to five times. 

To choose the number of repetitions, i.e. the duration of observations, for our targets we used the Standard Thermal Model (see 
SectionO to predict their flux densities in the PACS bands. Based on earlier Spitzer work (Stansberry et al.j 2008) we adopted a 



geometric albedo of 0.08 and a beaming parameter of 1.25 for observation planning purposes. The predicted thermal fluxes depend 
on the sizes, which are connected to the assumed geometric albedo and the absolute V-magnitudes via Equation ([3]). In some cases 
the absolute magnitudes used for planning purposes are quite different (by up to 0.8 mag) from those used for modeling our data as 
more recent and accurate visible photometry was taken into account (see Section l23] l. 

The PACS scan-map AOR allows the selection of either the blue or green channel; the red channel data are taken simultaneously 
whichever of those is chosen. The sensitivity of the blue channel is usually limited by instrumental noise, while the red channel 
is confusion-noise limited (PACS AOT release note 2010). The sensitivity in the green channel can be dominated by either source, 
depending on the depth and the region of the sky of the observation. For a given channel selection (blue or green) we grouped pairs 
of AORs, with scan orientations of 70° and 110° with respect to the detector array, in order to make optimal use of the rectangular 
shape of the detector Thus, during a single visit of a target we grouped 4 AORs to be observed in sequence: two AORs in different 
scan directions and this repeated for the second channel selection. 

The timing of the observations, i.e. the selection of the visibility window, has been optimized to utilize the lowest far-infrared 
confusion noise circumstances (IKiss et al. 20051 1 such that the estimated signal-to-noise ratio due to confusion noise has its maxi- 
mum in the green channel. Each target was visited twice with similar AORs repeated in both visits for the purpose of background 
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Table 2. Individual observations of the sample of 19 TNOs by Herschel/PACS. OBSIDs are the observation identifiers, duration is 
the total duration of the four AORs (see text), mid-time is the mean UT time, r is the mean heliocentric distance, A is the mean 
Herschel-target distance, and a is the mean phase angle. Each line corresponds to one visit and all three channels were observed at 
each visit. 



Target 


OBSIDs 




Duration 




Mid-time 




r 




A 


a 








(min) 












(AU) 


(AU) 


n 


119951 (2002 KX14) 


1342205144- 


■5147 


59, 


.8 


26-Sep^ 


■2010 


21 


:52: 


:54 


39, 


.3993 


39 


.8605 


1.30 


119951 (2002 KX14) 


1342205175- 


■5178 


59, 


.8 


27-Sep^ 


■2010 


15: 


:35: 


:54 


39, 


.3993 


39, 


.8716 


1.29 


(2001 XR„4) 


1342205184- 


■5187 


78, 


.6 


27-Sep^ 


■2010 


18: 


:48: 


:10 


44, 


.1004 


44, 


.4032 


1.25 


(2001 XR254) 


1342205264- 


■5267 


78, 


.6 


28-Sep^ 


■2010 


15: 


:48: 


:21 


44, 


.1004 


44, 


.3886 


1.25 


(2001 QY2,7) 


1342209492- 


■9495 


97, 


.4 


18-Nov^ 


■2010 


10: 


:06: 


:17 


43, 


.2452 


43, 


.3599 


1.31 


(2001 QY297) 


1342209650- 


■9653 


97, 


.4 


19-Nov^ 


■2010 


20: 


:50: 


:15 


43. 


.2455 


43, 


.3850 


1.31 


(2001 RZ141) 


1342199503- 


■9506 


97, 


.4 


01-Iuf 


■2010 


00: 


:52: 


:31 


41. 


.3006 


41, 


.6947 


1.31 


(2001 RZi4,) 


1342199614- 


■9617 


97, 


.4 


01-Iuf 


■2010 


20: 


:28: 


:08 


41. 


.3005 


41, 


.6819 


1.32 


(2002 GV^i) 


1342198847- 


■8850 


59, 


.8 


20-Iun^ 


■2010 


20: 


:07: 


:56 


40. 


.2818 


40, 


.5238 


1.42 


(2002 GV^i) 


1342198897- 


■8900 


59, 


.8 


21-Iun^ 


■2010 


21: 


:55: 


:30 


40. 


.2817 


40, 


.5412 


1.41 


79360 Sila" 


1342187073 




94, 


.4 


18-Nov^ 


■2009 


14: 


:24: 


:02 


43. 


.5090 


43, 


.2410 


1.27 


79360 Sila 


1342196137- 


■6140 


59, 


.8 


09-May^ 


-2010 


01 


:29: 


:23 


43. 


.5057 


43, 


.6530 


1.33 


79360 Sila 


1342196137- 


■6140 


59, 


.8 


10-May^ 


-2010 


09: 


:15: 


:12 


43. 


.5057 


43, 


.6753 


1.32 


88611 Teharonhiawako 


1342196099- 


■6102 


78, 


.6 


09-May^ 


-2010 


22: 


:19: 


:55 


45. 


.0856 


45, 


.3612 


1.24 


88611 Teharonhiawako 


1342196145- 


■6148 


78, 


.6 


10-May^ 


-2010 


18: 


:30: 


:58 


45. 


.0857 


45, 


.3475 


1.25 


120181 (2003 UR292) 


1342199618- 


■9621 


59, 


.8 


01-Iuf 


-2010 


21: 


:49: 


:00 


26. 


.7872 


27, 


.1477 


2.04 


120181 (2003 UR292) 


1342199646- 


■9649 


59, 


.8 


02-Iuf 


-2010 


12: 


:43: 


:03 


26. 


.7872 


27, 


.1379 


2.05 


(2005 EF298) 


1342208962- 


■8965 


97, 


.4 


03-Nov^ 


-2010 


20: 


:07: 


:09 


40. 


.7207 


41 


.1758 


1.25 


(2005 EF298) 


1342208999- 


■9002 


97, 


.4 


04-Nov^ 


-2010 


10: 


:08: 


:12 


40. 


.7207 


41 


.1667 


1.25 


138537 (2000 OKgy) 


1342197665- 


■7668 


78, 


.6 


03-Jun^ 


-2010 


00: 


:55: 


:51 


40. 


.3008 


40, 


.3581 


1.46 


138537 (2000 OK^) 


1342197717- 


■7720 


78, 


.6 


04-Iun^ 


-2010 


13: 


:55: 


:32 


40. 


.3006 


40, 


.3319 


1.46 


148780 Altjira 


1342190917- 


■0920 


76, 


.0 


21-Feb^ 


-2010 


23: 


:32: 


:17 


45. 


.5387 


45, 


.5571 


1.25 


148780 Altjira 


1342191120- 


■1123 


76, 


.0 


24-Feb^ 


-2010 


01: 


:31 


:06 


45. 


.5390 


45, 


.5940 


1.25 


(2002 KW14) 


1342204196- 


■4199 


59, 


.8 


09-Sep^ 


■2010 


11 


:11: 


:36 


40. 


.8385 


41, 


.0710 


1.38 


(2002 KW,4) 


1342204282- 


■4285 


59, 


.8 


10-Sep^ 


■2010 


09: 


:37: 


:01 


40. 


.8389 


41, 


.0870 


1.38 


(2001 KA77) 


1342205962- 


■5965 


78, 


.6 


06-Oct^ 


-2010 


21: 


:09: 


:43 


48. 


.1789 


48, 


.6245 


1.07 


(2001 KA77) 


1342206013- 


■6016 


78, 


.6 


07-Oct^ 


-2010 


14: 


:32: 


:45 


48. 


.1787 


48, 


.6354 


1.07 


19521 Chaos 


1342202285- 


■2288 


41, 


.0 


08-Aug^ 


-2010 


19: 


:40: 


:03 


41. 


.6914 


42, 


.1566 


1.24 


19521 Chaos 


1342202316- 


■2319 


41, 


.0 


09-Aug^ 


-2010 


11: 


:34: 


:49 


41. 


.6913 


42, 


.1464 


1.25 


78799 (2002 XWg^ 


1342190913- 


■0916 


57, 


.8 


21-Feb^ 


-2010 


22: 


:20: 


:43 


44. 


.6271 


44, 


.3163 


1.22 


78799 (2002 XW„) 


1342191116- 


■1119 


57, 


.8 


24-Feb^ 


-2010 


00: 


:19: 


:32 


44. 


.6279 


44, 


.3518 


1.23 


(2002 MS4) 


1342204140- 


■4143 


41, 


.0 


08-Sep^ 


-2010 


21: 


:12: 


:31 


47. 


.1673 


46, 


.9025 


1.19 


(2002 MS4) 


1342204292- 


■4295 


41, 


.0 


10-Sep^ 


-2010 


15: 


:03: 


:37 


47. 


.1670 


46, 


.9307 


1.20 


145452 (2005 RN43) 


1342195583- 


■5586 


41, 


.0 


25-Apr 


-2010 


22: 


:21 


:44 


40. 


,6885 


41, 


.1522 


1.26 


145452 (2005 RN43) 


1342195600- 


■5603 


41, 


.0 


26- Apr 


-2010 


21: 


:43: 


:12 


40. 


.6885 


41, 


.1378 


1.28 


90568 (2004 GV,) 


1342202869- 


■2872 


41, 


.0 


U-Aug^ 


-2010 


18: 


:35: 


:31 


39. 


.1876 


39, 


.3841 


1.46 


90568 (2004 GV,) 


1342202921- 


■2924 


41, 


.0 


12-Aug^ 


-2010 


15: 


:22: 


:58 


39. 


.1877 


39, 


.3984 


1.46 


120347 Salacia 


1342199133- 


■9136 


41, 


.0 


22-Iun^ 


-2010 


01: 


:17: 


:20 


44. 


.1464 


44, 


.0169 


1.33 


120347 Salacia 


1342199133- 


■9136 


41, 


.0 


22-Iun^ 


-2010 


18: 


:57: 


:52 


44. 


.1466 


44, 


.0058 


1.32 



Notes, r, A and a are from the JPL Horizons Ephemeris System ( [Giorgini et al. 1996] ). Observation in the chopped and nodded point-source 
mode and only at the blue and red channels JMuller et al.l 2010i. 



subtraction. The timing of the second visit was calculated such that the target has moved 30-50" between the visits so that the 
target position during the second visit is within the high-coverage area of the map from the first visit. Thus, we can determine the 
background for the two source positions. 

The observational details are listed in Table |2] All of the targets observed had predicted astrometric 3 cr uncertainties less than 
10" at the time of the Herschel observations (David Trilling, priv. comm.). 

2.3. Data reduction 

The data reduction from level (raw data) to level 2 (maps) was done using Herschel Interactive Processing Environment (HIP^l) 
with modified scan-map pipeline scripts optimized for the "TNOs are Cool" key program. The individual maps (see Fig. [T]for 
examples) from the same epoch and channel are mosaicked, and background-matching and source-stacking techniques are applied. 
The two visits are combined (Fig.|2]i, in each of the three bands, and the target with a known apparent motion is located at the center 
region of these maps. The detector pixel sizes are 3.2" x 3.2" in the blue and green channels, and 6.4" x 6.4" in the red channel 



' Data presented in this paper were analysed using "HIPE", a joint development by the Herschel Science Ground Segment 
Consortium, consisting of ESA, the NASA Herschel Science Center, and the HIFI, PACS and SPIRE consortia members, see 
|http : //herschel . esac . esa . int /DpHipeContribut ors . shtml , 
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Fig. 1. Individual maps of 120347 Salacia. Each map is the product of one observation (AOR). The first row (M1-M6) is from the 
first visit and the second row (M7-M12) from the follow-on visit. The first two columns (M1-M2, M7-M8) are observations in the 
100 fim or "green" channel and the others in the 160 i-im or "red" channel. The two scan angles are 1 10°(odd-numbered maps) and 
70°. The source is clearly seen in the map center in the green channel whereas the red channel is more affected by background 
sources and confusion noise. Orientation: north is up and east is to the left. 



whereas the pixel sizes in the maps produced by this data reduction are 1.1"/ 1 .4" / 2. 1 " in the blue / green / red maps, respectively. 
A detailed description of the data reduction in the key program is given in Kiss et al. (in prep.). 

Once the target is identified we measure the flux densities at the photocenter position using DAOPHOT routines (Stetson 1987]) 
for aperture photometry. We make a correction for the encircled energy fraction of a point source (P ACS photometer PSF 2010} for 
each aperture used. We try to choose the optimum aperture radius in the plateau of stability of the growth-curves, which is typically 
1.0-1.25 times the full-width-half-maximumof the PSF (5.2" / 7.7" / 12.0" in the blue / green / red bands, respectively). The median 
aperture radius for targets in the "TNOs are Cool" program is 5 pixels in the final maps (pixel sizes 1 . T'/l .4"/2. 1" in the blue / green 
/ red maps). For the uncertainty estimation of the flux density we implant 200 artificial sources in the map in a region close to the 
source (< 50") excluding the target itself. A detailed description of how aperture photometry is implemented in our program is 
given in lSantos-Sanz et an (2012). 

In order to obtain monochromatic flux density values of targets having a spectral energy distribution different from the default 
one color corrections are needed. In the photometric system of the PACS instrument ^mx density is defined to be the flux density that 
a source with a flat spectrum (AF^ ^constant, where A is the wavelength and is the monochromatic flux) would have at the PACS 
reference wavelengths ( Poghtsch et al. 20"T0) l. Instead of the flat default spectrum we use a cool black body distribution to calculate 
correction coefficients for each PACS band. The filter transmission and bolometer response curves needed for this calculation are 
available from HIPE, and we take as black body temperature the disk averaged day-side temperature calculated iterative ly for each 
target (| x Tss, using STM assumptions from Section[3] the Lambertian emission model and the sub-solar temperature from Eq.|2]) 
This calculation yields on the average 0.982 / 0.986 / 1.011 (flux densities are divided by color correction factors) for the blue / 
green / red channels, respectively, with small variation among the targets of our sample. 

The absolute flux density calibration of PACS is based on standard stars and large main belt asteroids and has the uncertainties 
of 3% / 3% / 5% for the blue / green / red bands (PACS photometer - Point Source Flux Calibration 201 1 1. We have taken these 
uncertainties into account in the PACS flux densities used in the modeling, although their contribution to the total uncertainty is 
small compared to the signal-to-noise ratio of our observations. 

The color corrected flux densities are given in Table |3] They were determined from the combined maps of two visits, in total 
4 AORs for the blue and green channels and 8 AORs for the red. The only exceptions are 19521 Chaos and 90568 (2004 GV9), 
whose one map was excluded from our analysis due to a problem in obtaining reliable photometry from those observations. The 
uncertainties in Table [3] include the photometric 1 cr and absolute calibration 1 cr uncertainties. 17 targets were detected in at least 
one PACS channel. The upper limits are the 1 cr noise levels of the maps, including both the instrumental noise and residuals from 
the eliminated infrared background confusion noise. 79360 Sila has a flux density which is lower by a factor of three in the red 
channel than the one published by IMiiller et al.l (2010). We have re-analyzed this earlier chopped/nodded observation using the 
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Table 3. Color corrected Herschel flux densities of the sample of 19 classical TNOs from coadded images of two visits. F70, Fioo 
and F160 are the monochromatic flux densities of the PACS blue / green / red channels. 



Target 




Fio 




-fioo 




-F160 




(mJy) 


(mJy) 




(mJy) 


119951 (2002 KX14) 


7.4 


±0.7 


10.2 


± 1.3 


7.2 


± 1.7 


(2001 XR2S4) 


2.5 


±0.7 




< 1.0 




< 1.4 


275809 (2001 QY297) 


1.1 


± 1.1 


4.2 


±0.8 


2.4 


± 1.1 


(2001 RZ143) 


2.2 


±0.7 




< 1.1 




< 1.1 


(2002 GV^i) 




< 0.8 




< 1.1 




< 1.7 


79360 Sila 


3.9 


±0.8 


5.5 


± 1.1 




< 5.7 


88611 Teharonhiawako 


1.9 


±0.7 


2.1 


±0.9 




< 2.1 


120181 (2003 UR292) 




< 4.6 


3.9 


± 1.1 




< 3.8 


(2005 EF298) 


1.4 


±0.6 


1.7 


±0.8 




< 1.9 


138537 (2000 OKgy) 




< 0.9 




<0.8 




< 3.1 


148780 Altjira 


3.4 


± 1.0 


4.3 


± 1.6 




< 2.1 


(2002 KW14) 


1.9 


±0.8 




< 3.0 




< 1.4 


(2001 KA77) 




< 2.1 


2.7 


± 1.0 




< 1.7 


19521 Chaos 


9.2 4 


:2.2° 


11.3 


± 1.2 


10.3 


± 1.8 


78799 (2002 XW93) 


17.0 


±0.9 


17.8 


± 1.5 


12.9 


± 1.9 


(2002 MS4) 


26.3 


± 1.3 


35.8 


± 1.5 


21.6 


±4.2 


145452 (2005 RN43) 


24.8 


± 1.3 


23.9 


± 1.8 


13.9 


± 1.9 


90568 (2004 GVg) 


16.9 


± 1.0 


19.3 4 


: 1.9* 


18.3 


±2.9 


120347 Salacia 


30.0 


± 1.2 


37.8 


±2.0 


28.1 


±2.7 



Notes. Observation 1342202317 excluded. <*> Observation 1342202923 excluded. 




Fig. 2. Combined maps of 120347 Salacia from the individual maps (Fig.[T]) in the green (left) and red channels. Orientation: north 
is up and east is to the left. 



latest knowledge on calibration and data reduction and found no significant change in the flux density values. As speculated in 
IMiiller et al.l (2010) the 2009 single-visit Herschel observation was most probably contaminated by a background source. 
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Table 4. Complementary Spitzer observations. The duration includes the total time of several visits. Observing epochs lasted 1-14 
days. The other quantities are as in Table|2] Targets below the horizontal line have / > 4.5°. 



Target 


Duration 


Starting day 


r 


A 


a 


Previous works 


This work 




(min) 


of observing 








MIPS F,4 


MIPS F70 


MIPS F04 


MIPS F70 






epoch 








(mJy) 


(mJy) 


(mJy) 


(mJy) 


119951 


128.10 


2005-08-26 


39.61 


39.59 


1.5 


0.0786 ± 0.0082" 


2.22 ± 1.45° 


0.080 ±0.013 


7.7 ± 1.7 


119951 


37.27 


2006-03-30 


39.58 


39.20 


14 


< 0.0363* 


< 3.90* 


0.083 ± 0.029 




275809 


594.87 


2008-11-21 


42.76 


42.39 


1.3 






<0.01 


< 1.5 


(2001 RZi4i) 


203.30 


2004-12-26 


41.38 


40.98 


1.3 


0.046 ± 0.008° 


< 0.7° 






79360 


494.70 


2008-05-18 


43.52 


43.30 


1.3 






0.05 ± 0.01 


< 2.7 


88611 


398.42 


2004-11-04 


45.00 


44.67 


1.2 






0.029 ± 0.010 


0.7 ± 0.6 


138537 


257.45 


2004-11-04 


40.57 


40.11 


1.3 


0.031 ±0.007" 


< 0.8° 






148780 


400.92 


2006-02-16 


45.33 


45.08 


1.2 


0.0167 ±0.0025° 


< 0.85° 


0.020 ± 0.008 


< 1.67 


(2002 KW14) 


213.05 


2005-08-26 


40.04 


40.03 


1.5 


< 0.006° 


3.3 ± 1.1° 






(2001 KA77) 


400.90 


2006-03-31 


48.56 


48.35 


1.2 


0.0077 ± 0.0023° 


4.12 ±0.81° 


< 0.025 


< 14 


19521 


66.06 


2004-09-24 


42.05 


41.68 


1.3 








9.0 ±3.0 


(2002 MS4) 


82.04 


2006-03-31 


47.39 


47.48 


1.2 


0.391 ±0.022* 


20.0 ±4.1* 


0.40 ± 0.02 


24.7 ± 2.9 


90568 


57.00 


2005-01-28 


38.99 


39.01 


1.5 


0.166 ±0.010* 


17.5 ± 2.2* 






120347 


227.62 


2006-12-03 


43.82 


43.39 


1.2 


0.546 ± 0.02F 


36.6 ± 3.7'- 







Notes. The Spitzer program IDs are: 55 (for the observations of Chaos), 3229 (Teharonhiawako), 3283 (2002 KX14 in 2006, 2002 MS4, 2004 
GV9), 3542 (2002 KX14 in 2005, 2001 RZ143, 2000 OK67, Altjira, 2002 KW14, 2001 KA77), 30 081 (Salacia), 4 0016 (Sila), and 50024 (2001 
QY297). References. *°' Color corrected MIPS flux density with 1 a measurement uncertainty from lBrucker et an (2009) together with quadrat- 
ically added calibration uncertainties. **' [Stansberry et al.|(2008), uncerta inty calculated from SNR and includes quadratically added calibration 
uncertainties, upper limits shown here are 1 cr. ^''' [Stansberry et al. (2012)| 



2.4. Complementary Spitzer observations 

About 75 TNOs and Centaurs in the "TNOs are Cool" program were also observed by the Spitzer Space Telescope 
jWerner et al. 2004l l using the Multiband Imaging Photometer for Spitzer (MIPS; IRieke et 1172 004). 43 targets were detected at 
a useful signal-to-noise ratio in both the 24 fim and 70 fim bands of that instrument. As was done for our Herschel program, many 
of the Spitzer observations utilized multiple AORs for a single target, with the visits timed to allow subtraction of background 
confusion. The MIPS 24 yum band, when combined with 70-160 jum data, can provide very strong constraints on the temperature of 
the warmest regions of a TNO. 

The absolute calibration, photometric methods and color corrections for the MIPS data are described in [Gordon et al. (200"7)1 



Engelbracht et al. (2007) and Stansberry et al. (2007) Nominal calibration uncertaintes are 2% and 4% in the 24 pm and 70 pm 
bands respectively. To allow for additional uncertainties that may be caused by the sky-subtraction process, application of color 
corrections, and the faintness of TNOs relative to the MIPS stellar calibrators, we adopt uncertainties of 3% and 6% as has been 
done previously for MIPS TNO data (e.g. Stansberry et al.| 2008, iBrucker et al. 12009). The effective monochromatic wavelengths of 
the two MIPS bands we use are 23.68 pm and 71.42 pm. With an aperture of 0.85 m the telescope-limited spatial resolution is 6" 
and 18" in the two bands. 

Spitzer flux densities of 13 targets overlapping our classical TNO sample are given in Table |4] The new and re-analyzed flux 
densities are based on re-reduction of the data using updated ephemeris positions. They sometimes differ by 10" or more from those 
used to point Spitzer. The ephemeris information is used in the reduction of the raw 70 pva data, for generating the sky background 
images, and for accurate placement of photometric apertures. This is especially important for the Classical TNOs which are among 
the faintest objects observed by Spitzer. Results for four targets are previously unpublished: 275809 (2001 QY297), 79360 Sila, 
88611 Teharonhiawako, and 19521 Chaos. As a result of the re processing of th e data, the fluxes for 11995 1 (2002 KX14), 148780 
Altjira, 2001 KA77 and 2002 MS4 differ from those pubUshed in lBrucker et"ari (2009) and [Stansberry et al.|(2008). 



The 70 pm bands of PACS and MIPS are overlapping and the flux density values agree typically within ±5% for the six targets 
observed by both instruments with SNR > 2. 



2.5. Optical photometry 

The averages of V-band or R-band absolute magnitude values available in the literature are given in Table [1] The most reliable way 
of determining them is to observe a target at multiple phase angles and over a time span enough to determine lightcurve properties, 
but such complete data are available only for 1 19951 (2002 KX14) in our sample. For all other targets assumptions about the phase 
behavior have been made due to the lack of coverage in the range of phase angles. 

The lAU (H,G) magnitude system for photometric phase curve corrections (IBowell et al. 1989] l. which has been adopted in some 



references, is known to fail in the case of many TNOs (Belskaya et al. 2008). Due to the lack of a better system we prefer linear 



methods as a first approximation since TNOs have steep phase curves which do not deviate from a linear one in the limited phase 



angle range usually available for TNO observations. The opposition surge of very small phase angles {a < 0.2°, Belskaya et al. 



2008) has not been observed for any of our 19 targets due to the lack of observations at such small phase angles. All of our Herschel 
and Spitzer observations are limited to the range 1 .0° < or < 2.1°. 
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Table 5. Overview of optical auxiliary data. The average absolute V-band (or R-band) magnitudes from literature are given in 
Table [T] and this table gives the absolute V-magnitudes with uncertainties which take into account the lightcurve (L.c.) amplitude 
AnjR (either in the UBVRI system R-band or Sloan's r'-band). The 0°-phase method (see text) tells how the extrapolation to zero 
phase was done (either in the reference or in this work). The last column indicates which corrections are significant in contributing 
to the corrected H^- 



Target 


Hv 


0° -phase method 


L.c. AniR 


L.c. period 


CoiTected Hy 


Corrections 




ref. 






(mag) 


(h) 


(mag) 


included 


119951 (2002 KX14) 


(a) 




phase study 






4.86 ±0.10 


assumed AniR 


(2001 XR254)* 


(b) 




(I) 


(> afew0.1:s)' 




6.05 ±0.15 


AmR 


275809 (2001 QY297)* 


(c,d,e) 




fl) 


0.49-0.66" 




5.86 ±0.31 


Aotr 


(2001 RZi43)* 


(f) 




R - ]fyf 






6.69 ±0.13 


assumed AniR 


(2002 GV31) 


(g) 










6.1 ±0.6 




79360 Sila* 


(h,i,j,k) 




(1) 


< 0.08" 




5.56 ± 0.04 


Aotr 


88611 Teharonhiawako* 


(1) 




(1) 


0.2" 


4.7526 ± 0.0007" 


6.00 ±0.13 


Aotr 


120181 (2003 UR292) 


(g) 




ySR = 0.14' 






7.4 ±0.4 




(2005 EF2,8) 


(g) 




(2) 






6.4 ±0.5 




138537 (2000 OY^-,) 


(m,n) 




= 0.15"70.14" 






6.47 ±0.13 


assumed AniR 


148780 Altjira* 


(o,p,d) 




(1) 


< 0.3" 




6.44 ±0.14 


AniR 


(2002 KW14) 


(q) 




(1) 


0.21/0.26^^ 




5.88 ±0.11 


Aotr 


(2001 KA77) 


(0) 




(1) 






5.64 ±0.12 


assumed Aotr 


19521 Chaos 


(h,i,n,r) 




(1) 


<0.10" 




5.00 ± 0.06 


AmR 


78799 (2002 XW,,) 


(g) 




(2) 






5.4 ±0.7 




(2002 MS4) 


(g) 




(2) 






4.0 ±0.6 




145452 (2005 RN43) 


(q) 




(1) 


0.04 ±0.01" 


5.62/7.32^ 


3.89 ± 0.05 


Aotr 


90568 (2004 GV,) 


(s) 


/J 


= 0.18 ±0.06'' 


0.16 ±0.03" 


5.86 ±0.03" 


4.23 ±0.10 




120347 Salacia* 


(q) 




(1) 


0.03 ± 0.01' 


6.09/8.1"^ 


4.24 ± 0.04 


AmR 



Notes. * denotes a known binary system JNoU et air20 08). (1) V-band linear phase coefficient yS = 0.14 ± 0.03 I Sheppard and Jewitt 
2002, /? calculated from values therein). (2) R-band linear phase coefficient j8r = 0.12 (Average from |Belskaya et al. 
20 08). References. |Rabinowitz et al. (2007)1 |Grundy et al. (2009)| IThirouin et akl (2012). [Grundy et al. (201 1)| 

[Doressound iram et al. (20 07)| ^' ISan tos-Sanz et al.' (2009). lAU Minor Planet Center / List of Tran sne ptunian Objects at 
fhttp : //www . minorplanetc enter . net/ i au/li st s/TNOs . html accessed June 2011. ^''^ [Boehnhardt et al. (2001)] 

JBarucci et al. (2000) [Grundy et a l. (2012) ~ " ' ' ' 

<°' Iporessoundiram et al. ( 2005)| ^"^"[^eppard (2007)] 
*'* Calculated from MPC data from 16 observations, 
component, amplitude calculated for the combined system. 



(0 



iBenecchi et al. (2009)[ 

Pema et al., in prep. 
"> [Sheppard and Jewittl (2002)! 
n. ^Dotto et al. (2008)] W|TE 



Benecchi et al.J (2011). 



jar 



Davies et al. (2000) 
[Doresso undiram et al. (2002) 



^gler and Romanishin (2000)[ DeMeo et al. (2009)_ 
[Osip et al. (2003)[ lightcurve caused by secondary 
irouin et al.l (2010). 



We have used the linear method (Hy = V - 5 log (rA) -fta, where V is the apparent V-magnitude in the Johnson-Coussins or the 
Bessel V-band, and other symbols are as in Table |2] and /3 is the lineaiity coefficient) to calculate the Hy values from individual V- 
magnitudes given in the references (see Table|5) for 2001 XR254, 275809 (2001 QY297), 79360 Sila, 88611 Teharonhiawako, 148780 
Altjira and 19521 Chaos. In order to be consistent with most of the Hy values in the literature we have adopted /3 = 0.14 ± 0.03 



calculated from Sheppard and Jewitt (2002). The effect of slightly different values of /3 or assumptions of its composite value used in 



previous works is usually negligible compared to uncertainties caused by lightcurve variability (an exception is 90568 (2004 GVg)). 

For five targets with no other sources available we take the absolute magnitudes in the R-band (//r in Table |5]) from the 
Minor Planet Center and calculate their standard deviation, since the number of V-band observations for these targets is very 
low, and use the average (V-R) color index for classical TNOs 0.59 ± 0.15 (Hainaut and Delsantil 2002) to derive Hy. While the 
MPC is mainly used for astrometry and the magnitudes are considered to be inaccurate by some works (e.g. Benecchi et al.i 201 1. 
[Romanishin and Tegler 2005), for the five targets we use the average of 9 to 16 R-band observations and adopt the average R-band 



phase coefficient /3r = 0.12 (calculated from Belskaya et al. 2008). 



The absolute V-magnitudes used as input in our analysis (the "Corrected Hy" column in Table |5]) take into account additional 
uncertainties from known or assumed variability in Hy. The amplitude, the period and the time of zero phase of the lightcurves 
of three hot classicals (145452 (2005 RN43), 90568 (2004 GV9) and 120347 Salacia) are available in the literature, but even for 
these three targets the uncertainty in the lightcurve period is too large to be used for exact phasing with Herschel observations. The 
lightcurve amplitude or amplitude limit is known for half of our targets and for them we add 88% of the half of the peak-to-peak 
amplitude (1 cr i.e. 68% of the values of a sinusoid are within this range) quadratically to the uncertainty of Hy. According to a study 
of a sample of 74 TNOs from various dynamical classes (Duffard et al. 200 9 ) 70% of TNOs have a peak-to-peak amplitude < 0.2 
mag, thus we quadratically add 0.09 mag to the uncertainty of Hy for those targets in our sample for which Ughtcurve information 
is not available. 
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3. Thermal modeling 

The combination of observations from thermal-infrared and optical wavelengths allows us to estimate various physical properties 
via thermal modeling. For a given temperature distribution the disk-integrated thermal emission F observed at wavelength A is 

F{A)^^-^ ^^B{A,T{S)) dS-M, (1) 

where e is the emissivity, A the observer-target distance, B{A, T) Planck's radiation law for black bodies, T {S) the temperature 
distribution on the surface S and u the unit directional vector toward the observer from the surface element dS. The temperature 
distribution of an airless body depends on physical parameters such as diameter, albedo, thermal inertia, and surface roughness. 

There are three basic types of models to predict the emission of an asteroid with a given size and albedo assuming an equi- 
librium between insolation and re-emitted thermal radiation: the Standard Thermal Model (STM), the Fast-rotating Isothermal 
Latitude thermal model jVeeder et al. 1989"] l. and the thermophysical models (starting from |Matson (1971)] e.g. Spencer et al. 1989] 



|Lagerros 1996 l. While originally developed for asteroids in the mid-IR wavelengths these models are applicable for TNOs, whose 
thermal peak is in the far-IR. 



The STM (cf. Lebofsky et al. 1986 and references therein) assumes a smooth, spherical asteroid, which is not rotating and/or 



has zero thermal inertia, and is observed at zero phase angle. The subsolar temperature Tss is 

\(1-A)Sq] 



Tss - 



(2) 



where A is the Bond albedo, 50 is the solar constant, rj is the beaming factor, cr is the Stefan-Boltzmann constant and r is the 
heliocentric distance. In the STM e does not depend on wavelength. The beaming factor rj adjusts the subsolar temperature. The 
canonical value i] = 0.756 is based on calibrations using the largest few main belt asteroids. The STM assumes an average linear 
infrared phase coefficient of 0.01 mag/degree based on observations of main belt asteroids. 

In this work we use the Near-Earth Asteroid Thermal Model NEATM ( Harris l9^8l l. The difference to the STM is that in the 
NEATM 77 is fitted with the data instead of using a single canonical value. For rough surfaces rj takes into account the fact that 
points on the surface S radiate their heat preferentially in the sunward direction. High values (77 > 1) lead to a reduction of the 
model surface temperature, mimicking the effect of high thermal inertia, whereas lower values are a result of surface roughness. 
Furthermore, the phase angle a is taken into account by calculating the thermal flux an observer would detect from the illuminated 
part of S assuming a Lambertian emission model and no emission from the non-illuminated side. 

Whenever data quality permits we treat 77 as a free fitting parameter However, in some cases of poor data quality this method 
leads to 77 values which are too high or too low and therefore unphysical. In these cases we fix it to a canonical value of 77 = 1.20+0.35 
derived by 'Stansberry et al. (2008) from Spitzer observations of TNOs. The physical range of 77 values is determined by using 
NEATM as explained in Mommert et aL(2012) to be 0.6 < 77 < 2.6. 

Throughout this work we assume the surface emissivity e - 0.9, which is based on laboratory measurements of silicate powder 
up to a wavelength of 22 yum (iHovis and Callahan 19661 ) and a usual approximation for small bodies in the Solar System. At far- 



IR wavelengths the emissivity of asteroids may be decreasing as a function of wavelength (Miiller and Lagerros 1998 1 from 24 to 
160 fim, but the amount depends on individual target properties. Thus, a constant value is assumed for simplicity. The variation 
of emissivity of icy surfaces as a function of wavelength could in principle provide hints about surface composition. H2O ice 
has an emissivity close to one with small variations jSchmitt et al. 1998] l whereas other ices may show a stronger wavelength 
dependence (e.g. [Stansberry et al. 1996j l. Near-IR spectroscopic surface studies have been done for five of our targets but none of 
them show a reliable detection of ices, even though many dynamically hot classicals are known to have ice signatures in their spectra 
jBarucci et al .'2011). 

For the Bond albedo we assume that A Ay since the V-band is close to the peak of the solar spectrum. The underlying 
assumption is that the Bond albedo is not strongly varying across the relevant solar spectral range of reflected light. The Bond 
albedo is connected to the geometric albedo via Ay - pvl, where q is the phase integral and pv the geometric albedo in V-band. 
Instead of the canonical value of ^ = 0.39 jBowell et al. 19891 ) we have adopted q - 0.336/?v + 0.479 (IBruckeret al.l 2009). which 
they used for a Spitzer study of classical TNOs. From the definition of the absolute magnitude of asteroids we have: 

PvSp,,;^7Ta^xl0ii'"O-"'\ (3) 

where Spi-oj is the area projected toward the observer, a is the distance of one astronomical unit, otq is the apparent V-magnitude of 
the Sun, and Hy is the absolute V-magnitude of the asteroid. We use niQ - (-26.76 + 0.02) mag jBessell et al. 1998] Hayes 1985] l. 
An error of 0.02 mag in Hy means a relative error of 1 .8% in the product ([3]). 

f4S~ 

We find the free parameters pv, D - J — ^ and 77 in a weighted least-squares sense by minimizing the cost function 



2 1 V' [F (Ai) - Fi^odeliAi)]^ 

= - y 2 ' w 

V ^ cr'r 



where v is the number of degrees of freedom, is the number of data points in the far- infrared wavelengths, F (/I,) is the observed 
flux density at wavelength /I, with uncertainty cr,, and F^odeX is the modeled emission spectrum. When only upper flux density limits 
are available, they are treated as having zero flux density with a 1 cr uncertainty equal to the upper limit flux density uncertainty 
(0 ± cr). The cost function (HJi does not follow the statistical distribution since in a non-linear fit the residuals are not normally 
distributed even if the flux densities had normally distributed uncertainties. 
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3.1. Error estimates 

The error estimates of the geometric albedo, the diameter and the beaming parameter are determined by a Monte Carlo method 
described in ' Mueller et al. (2 01"Ty' We generate 500 sets of synthetic flux densities normally distributed around the observed flux 
densities with the same standard deviations as the observations. Similarly, a set of normally distributed Hy values is generated. For 
those targets whose >> 1 we do a rescaling of the errorbars of the flux densities before applying the Monte Carlo error estimation 
as described in'Santos -Sanz et an (2012) and illustrated in iMommert et al. 1 (2012). 

The NEATM model gives us the eff'ective diameter of a spherical target. 70-80% of TNOs are known to be MacLaurin spheroids 
with an axial ratio of 1.15 (" Duff^ard et al. 20091 IThirouin et al. 1 2010). When the projected surface has the shape of an ellipse instead 
of a circle, then this ellipse will emit more flux than the corresponding circular disk which has the same surface because the Sun is 
seen at higher elevations from a larger portion of the ellipsoid than on the sphere. Therefore, the NEATM diameters may be slightly 
overestimated. Based on studies of model accuracy (e.g. Harris 20061 1 we adopt uncertainties of 5% in the diameter estimates and 
10% in the pv estimates to account for systematic model errors when NEATM is applied at small phase angles. 

4. Results of individual targets 

In this Section we give the results of model fits using the NEATM to determine the area-equivalent diameters (see Section |3]l as well 
as geometric albedos and beaming factors. We note that our observations did not spatially resolve binary systems, we therefore find 
area-equivalent diameters of the entire system rather than component diameters; this will be further discussed in Section lJ!2] 

In cases where also SpitzerfMlPS data are available for a target we determine the free parameters for both PACS only and the 
combined data sets. The solutions are given in Table|6] A floating-77 solution is only adopted if itsx^ is not much greater than unity. 
The exact limit depends on the number of data points. For N = 5 this limit is ^ 1-7. For the PACS-only data set we may adopt 
floating-/; solutions only if there are no upper-limit data points. 138537 (2000 OKev) has only upper limits from PACS, therefore 
only the solution using combined data is shown. 2001 RZ143 has five data points, three of which are upper limits. The PACS flux 
density at 70 fim is approximately a factor of three higher than the MIPS upper limit (see Tables |3] and |4|i. For this target we adopt 
the model solution determined without the MIPS 70 /im channel. The best solution for each target is shown in Fig.|3] 

The radiometric diameters determined with data from the two instruments are on the average close to the corresponding results 
using PACS data alone, but there are some significant differences as well, most notably 275809 (2001 QY297) and 2002 KW14. The 
former has a PACS-only solution, which is within the error bars of the green and red channel data points, but above the PACS blue 
channel data point. When the two upper limits from MIPS are added in the analysis the model solution is at a lower flux level below 
the PACS green channel data point but compatible with the other data (see Fig.|3]l. 2002 KW14 has upper limits in the PACS green 
and red channels as well as in the MIPS 24 fim channel. Without this upper limit in the shortest wavelength the PACS-only solution 
is at higher flux levels in short wavelengths and gives a lower flux at long wavelengths. 

When choosing the preferred solution we are comparing two fits with different numbers of data points used, therefore in this 
comparison we calculate for the PACS-only solution using the same data points as for the combined solution taking into account 
the different observing geometries during Herschel and Spitzer observations. In all cases where MIPS data is available the solution 
based on the combined data from the two instruments is the preferred one. 

2002 GV31 is the only non-detection by both PACS and MIPS in our sample of 19 targets. The astrometric 3 cr uncertainty at 
the time of the PACS observation was < 4" (semimajor axis of the confidence ellipsoicfl), which is well within the high-coverage 
area of our maps. 

The error bars from the Monte Carlo error estimation method may sometimes be too optimistic compared to the accuracy of 
optical data and the model uncertainty of NEATM (see Section[3]l. We check that the uncertainty of geometric albedo is not better 
than the uncertainty implied by the optical constraint (Eq.[3]l due to uncertainties in //y and niQ. The lower py uncertainty of four 
targets is limited by this (see Table|6]l. 

5. Discussion 

Our new size and geometric albedo estimates improve the accuracy of previous estimates of almost all the targets with existing 
results or limits (Table |7]i. Two thirds of our targets have higher albedos (Table |6]l than the 0.08 used in the planning of these 
observations, which has lead to the lower than expected SNRs and several upper limit flux densities (Table|3]l. Previous results from 
Spitzer are generally compatible with our new estimates. However, two targets are significantly different: 1 19951 (2002 KX14) and 
2001 KA77, whose solutions are compatible with the optical constraint (Eq.|3ll, but the new estimate has a large diameter and low 
albedo (target 1 19951) instead of a small diameter and high albedo, or vice versa (target 2001 KA77). It can be noted that there is a 
significant difference in the re-processed Spitzer flux densities at 70 ;um compared to the previously published values (see Table |4}, 
which together with the PACS 100 ;um data can explain the significant change in the diameter and geometric albedo estimates. Our 
new estimates for 148780 Altjira differ from the previous upper and lower limits based on Spitzer data alone (Table|7|i. This change 
can be explained by the addition of the 70 /vm and 100 /im PACS data points (Table O to the earlier MIPS 24 /vm value and the 
MIPS 70 yum upper flux density hmit (Tablegjl. 

The diameter estimates in our sample are ranging from 100 km of 120181 (2003 UR292) up to 930 km of 2002 MS4, which is 
larger than previously estimated for it. 2002 MS4 and 120347 Salacia are among the ten largest TNOs with sizes similar to those 
of 50000 Quaoar and 90482 Orcus. The size distribution of hot classicals in our sample is wider than that of the cold classicals, 

^ Asteroids Dynamic Site by A. Milani, Z. Knezevic, O. Arratia et al., [http : / /haitii 1 ton . dm . unipi . it/ astdys/ 1 accessed August 2011, 
calculations based on the OrbFit software. 
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Table 6. Solutions for radiometric diameters and geometric albedos (see text for explanations). For binary systems D is the area- 
equivalent system diameter In case of two solutions the preferred one is based on data from both PACS and MIPS instruments. 



Target 


Instruments 


D 

(km) 


Bin- 
ary? 


Pv 


U 


119951 (2002 KX14) 


PACS 

PACS, MIPS 


485!»^ 
455 ± 27 


no 


0.086--5 

0.097+"«" 

-0.013 


2 07+"-^** 

^•"'-0.82 

1.79+"!« 

-0.15 




PACS 


200+^" 


yes 


17+"" 


1 2 + 15 rfixed^ 


275809 (2001 QY,,,) 


PACS 

PACS, MIPS 


278+*^ 

-59 


yes 


104+""''* 

-0.050 

0.20+9 ?5 

-0.11 


1.2 ±0.35 (fixed) 
1.2 ± 0.35 (fixed) 


(2001 RZ,43) 


PACS 

PACS, MIPS* 


140+^' 

-33 


yes 


13+" ' "^ 

"■'--'-0.05 

0.191+9 

-0.045 


1.2 ± 0.35 (fixed) 
0.75+9 ?3 

-0.19 


(2002 GV31) 


PACS 


< 130 


no 


> 0.22 


1.2 ± 0.35 (fixed) 


79360 Sila 


PACS 

PACS, MIPS 


333+J^ 
343 ± 42 


yes 


095 
090+9 

"•"-^"-0017 


1.2 ±0.35 (fixed) 
1 36+9-?i 

^•''"-0.19 


88611 Teharonhiawako 


PACS 

PACS, MIPS 


234i^ 
177+^« 

-44 


yes 


n 19Q+0.062 

"•'^'^-o.oss 
0.22+"'^ 

-0.08 


1.2 ±0.35 (fixed) 
0.86+»^7 

-0.29 


120181 (2003 UR292) 


PACS 


104!g 


no 


16+" ''' 


1.2 ±0.35 (fixed) 


(2005 EF298) 


PACS 


174+27 

Vr 


no 


16+"" 


1.2 ±0.35 (fixed) 


138537 (2000 OKg,) 


PACS, MIPS 


151!^| 


no 


oW^ 

"•^"-0.08 


1.2 ±0.35 (fixed) 


148780 Altjira 


PACS 

PACS, MIPS 


313!| 
257i» 


yes 


0.048!»;»^| 


1.2 ±0.35 (fixed) 
1.46 ± 0.41 


(2002 KW14) 


PACS 

PACS, MIPS 


1811^ 

-81 


no 


24+" '■' 

"■^^-0.07 

0.08+9 

"-0.05 


1.2 ±0.35 (fixed) 
2.50+9" 

■ -1.64 


(2001 KA77) 


PACS 

PACS, MIPS 


252+*^ 
310+™ 

-60 


no 


155+""'" 
0.099+9 

-0.056 


1.2 ±0.35 (fixed) 
2.52+»l« 

-0.83 


19521 Chaos 


PACS 

PACS, MIPS 


600+'!!! 

-150 

600+''"' 

"""-130 


no 


050+" °'!° 

-^"-0.019 

050+'"'3"' 

"•"•'"-0016 


2 2+| -^ 
-1.0 

''•■^-11 


78799 (2002 XW53) 


PACS 


565+:?' 


no 


0.038+"°« 

-0.025 


0.79+"" 

-U.24 


(2002 MS4) 


PACS 

PACS, MIPS 


988!?^, 
934 ± 47'^ 


no 


046+""*^ 

051 +" 

"•"^^-0022 


1.2 ±0.35 (fixed) 
1.06 ± 0.06 


145452 (2005 RN43) 


PACS 


679^^^ 


no 


107+""^' 

"•'"'-0.018 


1.2 ±0.35 (fixed) 


90568 (2004 GV,) 


PACS 

PACS, MIPS 


670+2'" 

■-""-180 

680 ± 34^^ 


no 


0.08+55" 

0770+""™'' 

"■"' '"-0 0077 


1 9+' 5 

'•^^-0.8 

1 q'j+0.09 

^•"-0.07 


120347 Salacia 


PACS 

PACS, MIPS 


994+^^ 
901 + 45^ 


yes 


0361 +"""''' 

"•"-'"'-0.0055 

0.0439 ± 0.0044'' 


1 47+0.29 
'•^'-0.25 

1.156 ±0.031 



Notes. Lower uncertainty limited by the uncertainty of Hy for 275809 / PACS-MIPS, 2005 EF298, 78799, and 2002 MS4 / both solutions. 
(*) MIPS 70 fim channel excluded. ''^' Error estimate limited by the adopted diameter uncertainty of 5% of the NEATM model. '''' Error estimate 
limited by the adopted py uncertainty of 10% of the NEATM model. 

which are limited to diameters of 100-350 km (Fig.|4]i. The diameters of eight hot classicals from literature data (Table|7]) are within 
the same size range as the hot classicals in our sample. The cumulative size distribution of this extended set of 20 hot classicals 
(Fig.|5]) shows two regimes of a power law distribution with a turning point between 500 and 700 km. The slope of the cumulative 
distribution A^(> D) cc D '' is ^ » 1.4 for the 100 < D < 600 km (N=ll) objects. There are not enough targets for a reliable 
slope determination in the D > 600 km regime. The size distribution is an important property in understanding the processes of 
planet formation. Several works have derived it from the LF using simplifying assumptions about common albedo and distance. 
IFraser et al.l (2010) reported a slope of the differential size distribution of 2.8 ±1.0 for a dynamically hot TNO population (38 
AU < heliocentric distance < 55 AU and / > 5°). Our q + I based on a small sample of measured diameters of intermediate-size 
hot classicals is compatible with this literature value. The high-^ tail at D > 650 km in Fig. |5] indicates a change of slope when 
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Fig. 3. Adopted model solutions from Table |6] The black data points are from PACS (70, 100 and 160 fim) and the gray points are 
from MIPS (24 and 71 yum) normalized to the geometry of Herschel observations by calculating the NEATM solution (for given D, 
py and if) at both the epochs of the Herschel and Spitzer observations and using their ratio as a correction factor 



the population transitions from a primordial one to a collisionally relaxed population. Based on LF estimates, this change in slope 
for the whole TNO population was expected at 200 - 300 km (Kenyon et al. 2008 based on data from Bern stein et al. 2 004) or at 
somewhat larger diameters (IPetit et al 2006 1. For Plutinos a change to a steeper slope occurs at 450 km ( Mommert et al.l 2012). 

The dynamically cold and hot sub-populations are showing different geometric albedo distributions (Fig.|6) with the dynamically 
cold objects having higher geometric albedos in a narrower distribution. The average geometric albedo of the six cold classicals is 
0.17 + 0.04 (un-weighted average and standard deviation). The highest-albedo object is 8861 1 Teharonhiawako with pv = 0.22, or 
possibly 2002 GV31 with the lower limit of 0.22. These findings are compatible with the conclusions of Brucker et al.l (2009) based 
on Spitzer data that cold classicals have a high albedo, although we do not confirm their extreme geometric albedo of 0.6 for 1 1995 1 
(2002 KX14). 
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Fig. 3. continued. 



The darkest object in our sample is the dynamically hot target 78799 (2002 XW93) with a geometric albedo of 0.038. The 
highest-albedo hot classicals are found in the low-/ part of the sub-sample (see Fig. |7]l: 138537 (2000 OKe?) at / = 4.9° has 
pv = 0.20 and the inner belt target 120181 (2003 UR292) has a geometric albedo of 0.16. The 12 hot classicals in our sample have, 
on the average, lower albedos than the cold ones: /?v = 0.09 ± 0.05. The average of the combined hot classical sub-population of 
this work and Hterature is 0.1 1 + 0.04 if 55636 (2002 TX300) is excluded. 

From the floating-77 solutions of eight targets with data from both instruments included in the fitted solution (see Table |6] and 
SectionlU we have the average 77 = 1 .47 + 0.43 (un- weighted). Most of these eight targets have 77 > 1 implying a noticeable amount 
of surface thermal inertia. It should be noted, though, that inferences about surface roughness and thermal conductivity would 
require more accurate knowledge of the spin axis orientation and spin period of these targets. Our average t] is consistent with our 
default value of 1.20 ± 0.35 for fixed 77 fits. For comparison with other dynamical classes, the average beaming parameter of seven 
Plutinos is 77 = l.ll^g ( Mommert et al.l 2012) and of seven scattered and detached objects 77 - 1.14 + 0.15 CSantos-Sanz et al.1 
2012). The difference of using a fixed-77 = 1.47 instead of fixed-77 = 1.20 is that diameters would increase, on the average, by 10% 
and geometric albedos decrease by 16%. These changes are within the average relative uncertainties (19% in diameter and 57% in 
geometric albedo) of the fixed-?? solutions. 



5.1. Correlations 

We ran a Spearman rank correlation test ( [Spearman 1 904j to look for possible correlations between the geometric albedo py, diam- 
eter D, orbital elements (inclination /, eccentricity e, semimajor axis a, perihelion distance q), beaming parameter 77, visible spectral 
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Table 7. Adopted physical properties in comparison with previous works. Binary systems are indicated by 'B'. Column /Idetect lists 
the wavelengths used in the reference for radiometric diameters. If no radiometric result is available then binary system mass is used 
to give the diameter of the primary component assuming equal albedos. The second part of this table lists all other classical TNOs 
with radiometric or other reliable size estimates. The size of 50000 Quaoar is derived from both Spitzer and from direct imaging by 
Hubble. All targets in the lower part are dynamically hot. 



Target 



This work 
D (km) Pv 



n 097+0-014 

17+" 
20+''-^^ 

"■^"-0.11 
191+0066 

090+^-'^^^ 

"•"^"-0.017 

22+" ''' 
20+"^' 
0.07 1+or 



0.08!SJ;'^ 
0.099!""" 
050+''-^^^ 

"■"•^-0.016 

051+™ 
0770+" ""**■* 

'"-0.0077 

0.0439 ± 0.0044 



Previous work 
D (km) 



180+^'(| 
130-208" 
128-200" 

< 160 
250-420 
114-180" 

< 160 

< 200 
128-200" 

< 360 
634+^34 

< 742 
730 ± 120 

684+55 
720-1140" 
954 ± 80 



Pv 

60+"- * 

"•""-0.23 

0.09-0.23 
0.13-0.32 

> 0.23 
0.06-0.14 
0.13-0.32 

> 0.16 

> 0.10 
0.06-0.14 

> 0.05 

025+" "'" 

"•"^■^-0.008 

> 0.033 

073+" °' 

"•"'-'-0.03 

0.01-0.03 

0357+"-"°^^ 



Reference 



119951 (2002 KX14) 
(2001 XR2J4) 
275809 (2001 QY297) 
(2001 RZ143) 
79360 Sila 

88611 Teharonhiawako 
138537 (2000 OK67) 
148780 Altjira 

(2002 KW14) 
(2001 KA77) 
19521 Chaos 
(2002 MS4) 
90568 (2004 GV,) 
120347 Salacia 



455 ± 27 
200!*'^ 
200+1 
140+35 

343 ± 42 

1 '7'7+46 

i " _44 
151+3' 

257!»5 

3JO+170 
600+^^" 

"""-130 

934 ± 47 



680: 
901 : 



34 
45 



24, 71 
(binary) 
(binary) 
24 

70, 160 
(binary) 
24 
24 

(binary) 
71 

24, 71 
1200 
24, 71 
24, 71 
(binary) 
24, 71 



Brucker et al. (2009) 
Grundy etal. (2011) 



(2001 QD,,8) 
(1996 TS66) 
50000 Quaoar 
(2002 GJ32) 
20000 Varuna 
55637 (2002 UX25) 
55636 (2002 TX300) 
55565 (2002 AWi,,) 



24, 71 
24, 71 

24, 71, (direct) 

24, 71 

71 

24, 71 

(occultation) 
24, 71 



— 150^^" — 
190+50 

890 ± 70" 
220+^jJ 

680!!^" 
286 ± 10 
740 ± 100 



12+"-"' 

"•-^ -0.05 

0.18 ±0.04 
0.12!"-"' 
09+^^-8^ 

"•"^-0.03 

0.12t»'3 
0.88!":'^ 



Grun dy etal. (2011) 
Brucker et aL|(2009) 
MuUer et a lT|(2010) 
Grundy et al. (2011)| 
Brucker et aLj(2009) 
Brucker et al. (2009) 
Grun dy et al. (20TT)] 
Bru cker et al. 1 (2009) 
Bru cker et aLl(2009) 
Alten hoff etal. (2004) | 
Brucker etal. (2009) 
Brucker et al. (2009) 
Grundy et al. (2 011)] 
Stansb eiTy et al. (2012)| 
" Bruc ker et aLI (2009) 
Brucker et al.l(2009) 
Eraser and Brown (2010) 
LBrucker et al. (2009) 
Brucker et al. (2009) 
Stansberry et al.|(2 008) 
Ellio t et al. (2010)| 
Brucker et al.l(2009) 



Notes. Diameter of the primary component. 



Table 8. Selected correlation results (see text). 



Variables 


sub-sample 


Number of 
data points 


CoiTelation 
coefficient 


Significance 


Confidence 
limit (cr) 


DJ 


this work 


18 


69+" " 

"•"^-0.23 


0.002 


3.2 




this work / cold 


6 


-0.10+"''' 


0.8 


0.2 




this work / hot 


12 


n Q')+8:5S 
S-g+8:f? 


0.0011 


3.3 




this work and previous works 


26 


0.0004 


3.6 




this work / hot and previous works 


20 


60+" '* 

"•""-0 ?4 


0.005 


2.8 


D, Pv 


this work 


18 


-0 76+" " 
"• ' "-0.08 


0.0002 


3.7 




this work / cold 


6 


-0 40+"''** 

"•^"-0.41 


0.4 


0.8 




this work / hot 


12 


-0 62+"-2** 

"•"^-0.18 


0.03 


2.2 




this work and previous works 
this work / hot and previous works 


26 
20 


-0 59+""' 

-0 P 

-0 46+"-^' 


0.002 
0.04 


3.2 
2.1 


^v, 


this work 


18 


— 

"•-'"-o.is 


0.015 


2.4 




this work / cold 


6 


-0 02+"-''' 

"•" -0.58 


1.0 


0.03 




this work / hot 


12 


-0 74+<'-''5 

"• '^-0.18 


0.006 


2.7 


spectral slope, i 


this work / hot and previous works 


15 


-0 65+"-3'! 


0.011 


2.6 



slope, as well as B-V, V-R and V-I colorfl The Spearman correlation is a distribution-free test less sensitive to outliers than some 
other more common methods (e.g. Pearson correlation). We use a modified form of the test, which takes into account asymmetric 
error bars and corrects the significance for small numbers statistics. The details of our method are described in Peixinho et al. (2004)| 
and Santos-S anz et al.l (2012). The significance P of a correlation is the probability of getting a higher or equal correlation coef- 
ficient value -1 < p < 1 if no correlation existed on the parent population, from which we extracted the sample. Therefore, the 
smaller the P the more unlikely would be to observe a p if it was indeed equal to zero, i.e. the greater the confidence on the 
presence of a correlation is. The 99.7% confidence interval (3 cr), or better, corresponds to P = 0.003, or smaller. We consider a 
'strong correlation' to have |p| > 0.6, and a 'moderate correlation' to have 0.3 < |p| < 0.6. Selected results from our correlation 
analysis are presented in Tableland discussed in the following subsections. 



3 Erom the Minor Bodies in the Outer Solar System database, |http : / /www . eso ■ org/- ohainaut/MBOSS| accesses Nov 201 1. 
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Fig. 4. Distribution of diameters from this work (upper left), the cold classicals of this work (upper right), the hot classicals of this 
work (lower left), and all hot classicals including literature results from Table |7](lower right). The last plot includes only dynamically 
hot classicals. The bin size is 100 km. 
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Fig. 5. Cumulative size distribution of dynamically hot classicals from this work and literature (Table |7}. The power law has a 
change between 500 and 700 km. The intermediate size classicals have a slope parameter of q - 1 .4. 



5.1 .1 . Correlations with diameter 

We detect a strong size-inclination correlation in our target sample (see Fig. |7] and Table |8]l. When Uterature targets, all of whom 
are dynamically hot, are included in the analysis we get a correlation of similar strength. Previously this presumable trend has 
been extrapolated from the correlation between intrinsic brightness and inclination (iLevison and Steriil 2001). We see this strong 
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Fig. 6. Distribution of geometric albedos from this work (upper left), the cold classicals of this work (upper right), the hot classicals 
of this work (lower left), and all hot classicals including literature results from Table Q (lower right). The bin size is 0.05. The 
Haumea family member 55636 (2002 TX300) with py - 0.88 is beyond the horizontal scale. 



size-inclination correlation also among the hot classicals sub-sample, but not among our cold classicals where we are limited by the 
small sample size. 

Other orbital parameters do not correlate with size. We find no correlation between size and colors, or spectral slopes, nor 
between size and the beaming parameter 77. The possible correlation between size and geometric albedo is discussed in Section l5.1.2l 



5.1 .2. Correlations with geometric albedo 

We find evidence for an anti-correlation between diameter and geometric albedo, both in our sample and when combined with 
other published data of classical TNOs (see Fig. [8] and Table |8]l. Other dynamical populations with accurately measured diame- 
ters/albedos show a different behavior: there is no such correlation seen among the Plutinos jMommert et al.l 2012) and a combined 
sample of 15 scattered-disc and detached objects show a positive correlation between diameter and geometric albedo at 2.9 cr level 
aSantos-Sanz et al.l 2012). 

As it might be suggested visually by the distribution of diameters of classical TNOs (see Figs.|4]and[8]l, we have analyzed the 
possibility of having two groups with different size-albedo behaviors, separating in size at Z) a; 500 km regardless of their dynamical 
cold/hot membership. We have found no statistical evidence for it. 

With our method of accounting for error bars, which tend to 'degrade' the correlation values, geometric albedo does not correlate 
with Hv, orbital parameters, spectral slopes, colors, nor beaming parameters 77. Also when the literature targets are added we find 
no evidence of correlations. 



5.1.3. Other correlations 

The known correlations between surface color/spectral slope and orbital inclination ( Trujillo and Brown 2002 IHainaut and Delsantil 
2002), and between intrinsic brightness and inclination (ILevison and Stenil 2001). usually interpreted as a size-inclination correla- 
tion, might lead us to conclude there was a consequent color/slope-size correlation. Our analysis with measured diameters does 
not show a correlation neither with spectral slope nor visible colors, as one might expect. Note, however, that we do not possess 
information on the surface colors/slopes of ~ 1/2 of our targets (~ 1/3 when complemented with other published data) leading 
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Fig. 7. Radiometric diameter as well as geometric albedo vs inclination. The cold classicals of our sample are marked with blue 
squares, hot classicals with red crosses and hot classicals from literature with gray color The high-albedo target 55636 (2002 TX300) 
is beyond the scale. 



17 



E. Vilenius et al.: "TNOs are Cool": A survey of the trans-Neptunian region 



0.4- 




100 200 400 600 800 1000 

Radiometric diameter (km) 



Fig. 8. Geometric albedo vs radiometric diameter (blue triangles - cold classicals, red crosses - hot classicals from our sample, 
gray = other hot classicals from literature, see Table |7}. 

to the non-detection of the color/slope-inclination trend, which is known to exist among classicals. Thus, a more complete set of 
color/slope data would be required for our targets. Only when combining the hot sub-sample from this work and literature we see 
a non-significant anti-correlation between slope and inclination (see Table [8]l. We do not find any correlations of the B-V, V-R and 
V-I colors with other parameters. 

The apparent Hy vs / anti-correlation in our target sample mentioned in Section ITT] is almost significant (2.1 cr) for our hot 
sub-population (see Tablets). 



5.2. Binaries 

Binary systems are of particular scientific interest because they provide unique constraints on the elusive bulk composition, whereas 
all other observational constraints of the composition only pertain to the surface of the object. The sizes of binaries can be con- 
strained based on the relative brightness difference of the primary and the secondary components, but only if suitable assumptions 
about the relative geometric albedo are made. Alternatively, geometric albedos can be constrained under certain assumptions about 
the relative sizes. The ranges given in the literature are usually based on the following assumptions: i) the primary and secondary 
objects are spherical, ii) the primary and secondary have equal albedos, and iii) objects have densities within a limited assumed 
range. 

Six of our targets are binaries with known total mass m and brightness difference between the two components AV (see Tables- 
Assuming the two components to have identical albedos, the latter can be converted into an area ratio and, assuming spherical shape, 
a diameter ratio k - D2/D1 (with component diameters Di and D2). Component diameters follow from the measured NEATM 
diameter D (see Table|6]l and D2/D1: - D\ + D\ (since D is the area-equivalent system diameter). This leads to a "volumetric 

diameter" D' = ^^^=£>- The mass densities are given in Table|9] Within the uncertainties, the measured bulk densities scatter 
around roughly 1 g cm"^, consistent with a bulk composition dominated by water ice, as expected for objects in the outer Solar 
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Table 9. New density estimates for binaries. 



Target 


Adoptea AV 


Mass*^ 


Bulk density 




(mag) 


(xlO"*kg) 


(g cm-^) 


(2001 XR254) 


0.43 


4.055 ± 0.065 


1 4+'-^ 


275809 (2001 QY297) 


0.20 


4.105 ±0.038 


lJ|| 


79360 Sila 


0.12* 


10.84 ± 0.22 


0.73 ±0.28 


88611 Teharonhiawako 


0.70 


2.445 ± 0.032 


I 14+0.87 


148780 Altjira 


0.23 


3.986 ± 0.067 




120347 Salacia 


2.32 


466 ± 22 


1.38 ±0.27 



Notes. References. [Grundy et al. (201 1)| "" [Grundy et al. (2012)| 

System. Significant mass contributions from heavier materials, such as silicates, are not excluded however, and would have to be 
compensated by significant amounts of macroporosity. The largest object, 120347 Salacia, has a bulk density > 1 g cm"-*. This could 
indicate a lower amount of macroporosity for this object, which is subject to significantly larger gravitational self-compaction than 
our other binary targets. 

6. Conclusions 

The number of classical TNOs with both the size and the geometric albedo measured radiometrically is increased by eight from 
22 to 30. Four other targets, which previously had estimated size ranges from the analysis of binary systems, now have more 
accurate size estimates. The number of targets observed and analysed within the "TNOs are Cool" program (Miilleret al.'2010, 
ILellouch et al.l 2010. ILim et al.l 2010. ISantos-Sanz et al.l 2012. IMommert et al.l 2012) is increased by 18 and the observation of 
79360 Sila disturbed by a background source in Mii ller et al. 1 (20 10) has been re-observed and analyzed. Furthermore, three targets 
which earlier had upper and lower limits only based on Spitzer data alone (148780 Altjira, 138537 (2000 OKe?) and 2001 RZ143) 
now have accurately estimated diameters and albedos. The new Altjira solution is outside of the previous limits based on Spitzer 
data alone. The three PACS data points near the thermal peak are providing reliable diameter/albedo solutions, but in some cases 
adding Spitzer data, especially the 24 pm data point in the lower-wavelength regime, constrains the solution and allows smaller 
error bars and more reliable estimates of the beaming parameter Compared to previous works the size estimates of 119951 (2002 
KX14), and 2002 MS4 have increased. 2002 MS4 (934 km) is similar in size to 50000 Quaoar and the refined size of 120347 Salacia 
(901 km) is similar to that of 90482 Orcus. We find a diameter for 2001 KA77, which is approximately half of the previous estimate 
jBrucker et al. 2009), and a geometric albedo approximately 4 times higher The largest change in estimated geometric albedo is 
with 119951 (2002 KX14) from 0.60 to 0.097. 

The main conclusions based on accurately measured classical TNOs are: 

1. There is a large diversity of objects' diameters and geometric albedos among classical TNOs. 

2. The dynamically cold targets have higher (average 0.17 + 0.04) and differently distributed albedos than the dynamically hot 
targets (0.09 + 0.05) in our sample. When extended by seven hot classicals from literature the average is 0.1 1 ± 0.04. 

3. Diameters of classical TNOs strongly correlate with orbital inclination in the sample of targets, whose size and geometric albedo 
have been accurately measured, i.e. low inclination objects are smaller We find no clear evidence of an albedo-inclination trend. 

4. Our data suggests that geometric albedos of classical TNOs anti-correlate with diameter, i.e. smaller objects possess higher 
albedos. 

5. Our data does not show evidence for correlations between surface colors, or spectral slope, of classical TNOs and their diameters 
nor with their albedos. 

6. We are limited by the small sample size of radiometrically measured accurate diameters/albedos of dynamically cold classicals 
(N=6) finding no statistical evidences for any correlations. 

7. The cumulative size distribution of hot classicals based on the sample of measured sizes in the range of diameters between 100 
and 600 km (N= 1 1 ) has a slope of ^ a; 1 .4. 

8. We determine the bulk densities of six classicals. They scatter around ~ 1 g cm"-'. The high-mass object 120347 Salacia has a 
density of (1.38 ± 0.27) g cm^^ 
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